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Objectives 

We shall 

• review the maths needed, 

• derive Maxwell's equations, 

• derive the wave equation, 

• derive waveguide modes, 

• derive synchrotron radiation, 

• review vector potential, and 

• derive the Hamiltonian. 
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Wave Equation 



In a stretched string, tension is T and mass per unit length is p. 
Consider a small section dx that is displaced. 



T 




resultant vertical force = mass x vertical acceleration 

a, 



Ty {x + dx) — Ty (x) = pdx x 



dt 2 



Divide both sides by dx, 



T d 2 y = d 2 y 
dx 2 dt 2 



Solution has a very general form (can check by substution) 



y = f(x — vt) where v = 



T 
P 
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Wave Equation 



is a solution, where v = JT/p. This means a function of any 
shape can travel along the string - just what we expect when 
we shake one end of a string, v is then the velocity. 

For particular frequency uu = 2irf and small displacements, we 
often assume 

y = sin(/c(x — vt)) — s\r\(kx — out) where v = — 

k 

In terms of v, the equation is 

d 2 y _ 1 d 2 y 
dx 2 v 2 dt 2 
In 3D, the equation generalises to 

d 2 f | d 2 f | d 2 f = 1 d 2 f 

dx 2 dy 2 dz 2 v 2 dt 2 

where f(x,y,z) is some 3D function which could be pressure, 
electric field, etc. 
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Differential 



A(x,y,z,t) is a function of many variables. The change in A as 
a result of small changes in the variables is the differential 

a A = ax H ay -\ dz -\ at 

dx dy dz dt 

If we need the time rate of change of A when all variables are 
changing, we cannot just take dA/dt. We must divide all of dA 
by dt: 

dA dAdx . dAdy dAdz dAdt 

— = 1 -{ 1 

dt dx dt dy dt dz dt dt dt 

We often use x for dx/dt, so 

dA dA. , dA. , dA. , dA 

— — — x H y H 2; H 

(it 9x 9i/ dz dt 

Note that ^ and ^ are completely different. In x, y and z 
must be fixed. 
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Vector Calculus - Definitions 



Scalar function (or scalar field) 



</>(x,y,z) or 0(x) 
Vector function (or vector field) 

A(x, y, z) = iA x (x, y, z) + ]A y (x, y, z) + kA z (x, y, z) 
Gradient 

ax ai/ 



Curl 



V x A = 



i 

dx 

Ax 



J k 



dy dz 

Ay A z 



Divergence 



_ A dA x , 9A y 9A Z 
V.A = 1 H 

<9x dy dz 
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Vector Calculus - Identities 



Divergence of a curl is zero 

V.(V x A) = 0 
Curl of a gradient is zero 

V x (V0) = 0 

Curl of a curl 

V x (V x E) = V(V.E) - V 2 E 
Divergence of cross product 

V.(E x H) = (V x E).H - (V x H).E 

To find more, google "vector calculus identities". 

You can prove these using the definitions of div, grad and curl. 
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Divergence - Physical Meaning 



Relates the vector field in enclosed in a volume to the flux of 
the field at the surface. Useful for electric and magnetic fields. 

Start with physical meaning of divergence. 



A 




Define flux of a field A through a small area dS as 

flux = A.dS 

The dot product means the same as projecting dS in the 
direction of A so that 

flux = AdS cos 0 
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Divergence - Physical Meaning 



-A x (x, y, z) 




dz 



> A x (x l dx, y, z) 



dy 



Consider a small volume element dxdydz. The sum of flux at 
two opposite sides along x direction are 

dA x 

A x {x + dx, y, z)dydz — A x {x, y, z)dydz = dxdydz 

dx 

The total flux over all sides is therefore 

7C , dA x dA v dA z , , 

A.dS = dxdydz -\ dxdydz -\ dxdydz (1) 

SS dx dy ' dz 

This gives the physical meaning of divergence: 

^ = V.A (2) 

Divergence at a point is the net flux coming out of the point 
per unit volume. 
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Divergence Theorem 



Now imagine a large volume V enclosed by surface S. This is in 
a vector field A. Divide this volume into small volume elements 
packed together. 




Consider the flux of A from each element. Clearly, fluxes from 
touching surfaces cancel. 
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Divergence Theorem 



This means that all fluxes cancel except those at the outermost 
surface S. 

Recall the flux at each element from Eq. (1) 

A.dS = — — -dxdydz + — — -dxdydz + — — -dxdydz 
5S dx dy dz 

Summing over all elements gives 

A.dS = [ V.AdV 
This is the divergence theorem. 
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Curl - Physical Meaning 



Consider a loop C in a vector field A. Consider a short 
segment dl along C. 




At each segment, take the dot product A.dl. If we think of A 
as force, then A.dl is work done. Sum this up over all segments 
of the loop. This is called the 

circulation = / A.dl 

Jc 

If the loop becomes very small, the circulation approaches the 
curl component normal to the loop. 
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Curl - Physical Meaning 



To see this, project a small, square loop onto the yz plane 
through the loop. 



sc 




- y 



Contribution from sides z direction is 



A z (x, y + dy, z)dz — A z (x, y, z)dz = dydz 

dy 
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Curl - Physical Meaning 



Adding the contribution from sides along y direction gives 

— — -dydz — — — -dzdy = (V X A) x dydz 
dy dz 

which is the x component of curl A. 



Adding all 3 components together gives 

/ A.d\= (V x A).5S (3) 

JSC 
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Stoke's Theorem 



Now consider a larger loop C. Imagine a surface enclosed by 
this loop. It can be any of the possible surfaces that has C as 
boundary. 



Divide S into small square elements SS. 

Add up the circulations due to all elements. Contributions from 
adjacent elements to circulation along a common edge cancel. 
Only contributions along the outermost C remain. 

Summing Eq. (3) over all elements gives Stoke's theorem: 





14 



Gauss's Law 



Consider a volume V enclosed in surface S. We want to show 
that the total flux of electric field through S is related to the 
total charge in V. 




Suppose there is a point charge Q in V. Consider a narrow 
cone of field from Q crossing the surface at dS. 

Area dS varies as r 2 from Q. Field E varies as 1/r 2 So flux is 
the same for any surface element in the cone. 
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Gauss's Law 



Therefore the total flux is equal to that through a sphere 
centre at Q. 




So total flux through S is 

E.dS = Q 0 x 47rr 2 = ^ 

'S 4-7T€Qr Z €Q 

For arbitrary distribution of charges p(r), Q generalises to the 
total charge enclosed: 

1 



E.dS = 



p(r)dV 
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Gauss's Law 



Using divergence theorem: 

V.EdV = — [ o(r)dV 

'V 6Q JV 

Because this is true for any volume V, the integrands must be 
equal: 

This and the integral above are the two forms of Gauss's law. 
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Gauss's Law for Magnetism 



A magnet always has north and south pole. The field 
distribution from each pole is like that of an electric charge. So 
Gauss's law also applies to magnetic field. 

Because north and south poles are like opposite charges, net 
"magnetic" charge is always zero. Gauss's law for magnetism 
can be written down from Gauss's law for electricity. Just 
replace E by B and charges by zero: 




V.B = 0 
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Ampere's Law 



A current density I through a straight conductor produces a 
circular magnetic field B around it related by Ampere's law: 



/ B.dl = Linl 

Jc 



c 



Generalising to current density J, 

/ B.dl = Lin J.dS 

Jc Js 

This is true for a surface through a current. But what if S goes 
through the empty space of a capacity with no J? 
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Ampere's Law 




Maxwell generalised J to include changing E field between 
capacitor. E is related to charge Q on capacitor of area A by 

£ = ^(From Gauss's law. How?) 

e 0 A 

This is related to current by 

dQ A dE 

I = -f - = e Q A 

dt dt 

and to current density by 

I dE 

J = — = en 

A U dt 
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Ampere's Law 



This is called displacement current. Ampere's law is then 
generalised to 

( B.dl = an [ J.dS + enLio [ — .dS 

Jc pu Js upu Js at 

Applying Stoke's theorem, 



J V x B.dS = mo J J-dS + e 0 ^o J ^ 



dS 



Since this is true for and surface S, the integrands on two sides 
must be equal: 

V x B = i_i 0 J + cquo 



dt 

This and the above integral equation are the 2 forms of 
Ampere's law with displacement current. 
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Faraday's Law 



When the magnetic flux 0 linking a coil changes, a voltage V is 
induced round the coil C: 




N S 



In terms of magnetic field B and electric field E\ 

[ E.dl = -— [ B.dS 

JC dt JS 
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Faraday's Law 



Using Stoke's theorem, 

, V x E.dS = — [ B.dS 

'S dt JS 

Since this is true for any surface S, the integrands must be 
equal: 

dB 
V x E = — 



dt 

This is the differential form of Faraday's law, 
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Magnetic Force 



Consider a horizontal rectangular wire loop with one sliding side 
and vertical magnetic field through it. 





V 



A 



x 



A A A A A A A A A A A A A A A A B 



In terms of magnetic field B and electric field As the side 
moves with velocity v, flux changes and voltage V is induced 
Faraday's law leads to 



V = - 



d(BA) d(Blx) 



= Blv 



dt dt 

Current induced can be understood as due to forces on 
electrons in moving side. 
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Magnetic Force 



V is work done w per unit charge q. Work done is due to force 
F acting on charge over side /: 




q / F 



A A A A A A A A A 



w/q = Fl/q = Blv =>* F = qBv 

If B is at some angle 9 to v, only component of B 
perpendicular to v has effect. So 

F — qBv Sin 0 =^> F = qv x B 



Including the force from any electric field E, 

F = qE + qv X B 

This is resultant is called the Lorentz force. 
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Wave Equation 



Maxwell's equations: 



Gauss's law 



Gauss's law for magnetism 



Ampere's law 



Faraday's law 



V.B = 0 
V x B = fio J + cquo 



dE 

~dt 



V x E = - 



<9B 

~at 



In a vacuum with no charge and current, these become 



Gauss's law 

Gauss's law for magnetism 
Ampere's law 

Faraday's law 



V.E = 0 
V.B = 0 

V x B = e 0 MO 



dE 



dt 



V x E = - 



dB 
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Wave Equation 



Take the curl of Faraday's law: 

d 

V x (V x E) = V x B 

dt 

Apply a vector identity to the left and Ampere's law to the 
right: 

3 <9E 

V(V.E) - V 2 E = -eono-— 

at at 

Using Gauss's law, first term on left is zero. So 

2-r. 9 2 E 
V E = eoMo^ 

E.g. x component is: 

d 2 E x d 2 E x d 2 E x ^ d 2 E x 

dx 2 dy 2 dz 2 dt 2 
Recall the 3D wave equation 

d^l + d 2 i + d 2 i = 1 d 2 f 

dx 2 dy 2 dz 2 v 2 dt 2 
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Wave Equation 



Clearly, we had a wave equation for electric field. The solution 
could be any electric field variation that travels with speed v. 
Comparing the equations, 

v = 1 = 2.998 x 10 8 m/s 

the speed of light! 

1864: This result from Maxwell is strong evidence that light is 
an electromagnetic wave. 

1888: Hertz generated a new (at that time) type of wave 
which we now call microwave - by using spark gap to create 
disturbance in electric field. 

Any waveform travelling at speed of light is a solution. For 
wave of a single frequency, consider 

E x = E 0 COS(k x x - cut) 
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Wave Equation 



For 3D, we can check by substitution that this is also a solution 

E x = Eq COS(k x x + k y y + k z z — out) = Eq COS(k.r — uot) 

This is a wave travelling in the k= (k x ,k y ,k z ) direction. It is 
real part of the complex form: 

E = E 0 e i(k - r - wt) 

where the E y and E z are included to form E. Substituting into 
Faraday's law and doing the differention gives 

ik x Eq = ic^Bo 

1. So an electric wave must generate a magnetic wave. 

2. The magnetic oscillation must be perpendicular to the 
electric one. 

3. Both magnetic and electric oscillations must be 
perpendicular to the wave direction k. 

This describes the transverse electromagnetic wave in free 
space. 



29 



Waveguide 



Waveguides are used to transmit electromagnetic waves with 
minimal loss. A familiar form is the long, uniform, hollow 
metallic pipe used for microwave. 



Suppose the axis is along z. For a fixed frequency, field H 
contains factor e~ lujt . Time derivative gains factor — iu. So 
Faraday's law becomes: 



V x E = - 



dt 



• 

1 


• 

J 


k 


d 


d 


d 


dx 


dy 


dz 


E x 


Ey 


E z 



here /xH = B. In full 



dE z dE 



y _ 



dy 
dE z 

dx 
dE 



dz 
dE x 



X 



y 



dz 
dE 



X 



dx 



dy 
30 



iuoiiHy 



Waveguide 



We drop the subscripts to eo and /xq so the equations also apply 
to wave in a medium. 

Assume (as trial solution) that spatial variation of z is e lkzZ . So 

dE z 



— ik z E y = iu)fiH x 

dy 



dE 



z 



— ik z E x ^j — iuoiiH 



dx 

dEy dE x 
dx dy 



y 



A similar expansion of Ampere's law gives 

dH z 



ik z H y = —iuj€E x 

dy 



dH 



z 



dx 



— ik z H x ) — —iujeEy 



dHy dH x 
dx dy 



— —iujeE: 
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Waveguide 



Solving for transverse components: 



i ( 8E Z dH z 

tlx = -77^ cje— k z 



\ dy dx 



H 1 (ue— k? Hz 
y k£ \ dx dy 

i (, 8E Z . dH z 

Ex = T2\ kz ^ — l~ u ^^r 

k£ \ ox ay 
i ( dE z dH z 
k£ \ oy ox 

where k 2 = uj 2 hc — fcf . 

ALL TRANSVERSE COMPONENTS OF E and H CAN BE 
DETERMINED FROM ONLY E z and H z . 

Since E z and H z each have its own wave equation, we can solve 
one assuming the other is zero and vice versa. The separate 
solutions can then be added up to give a general solution. 
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Waveguide 



Then the solutions divide naturally into two families: 

1. When E z = 0 - Transverse Electric (TE) 

2. When B z = 0 - Transverse Magnetic (TM) 

To solve the equations, we need boundary conditions, assuming 
perfect conductor: 

1. E perpendicular to surface. So E z = 0. (Used to solve E z 
equation. So B z = 0 - TE.) 

2. B parallel to the surface. So dB z /dn = 0 - n means normal 
to surface, e.g x or y. (Used to solve B z equation. So E z = 0 - 
TE.) 
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Waveguide 



E.g. consider TE in rectangular waveguide. Wave equation for 
B z as same form as E x equation from earlier slide: 

d 2 B z , d 2 B z , d 2 B z d 2 B z 



<9x 2 ch/ 2 <9z 2 ' <9t 2 

Try separable variables: 

B z = B 0 (x i y)e ikz -" t 

Dropped subscript for k z as z is only propagating direction 
possible. Substituting into wave equation gives 

d 2 d 2 

" + tt^ + M^ 2 - fc I B z = 0 



dx 2 dy 2 
Applying boundary condition, find 



B z = H 0 cos — 

V a 



( m7rx\ ( Tiny 



COS 



b 

where m,n are any integers. Substituting into wave equation 
gives 

2 ,2 2 / m2 , n2 

lieuj — k — 7r — -+- 



a 



2 ' ft 2 
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Waveguide 



For s\r\(kz - out) to be a wave, \k\ > 0. When k = 0, this 
equation gives the smallest uj - cutoff frequency. 

Wave for each possible set of m, n is called a mode. 



E.g. a = 26, cutoff frequencies ratios (Jackson, 3rd ed.): 



\. n 
m \^ 


0 


1 


2 


3 


0 




2 


4 


6 


1 


1 


2.24 


4.13 




2 


2 


2.84 


4.48 




3 


3 


3.61 


5 




4 


4 


4.48 


5.66 




5 


5 


5.39 






6 


6 









In a waveguide, waves can be produced using a small antenna 
or coupler that sticks into the guide. Any mode with cutoff 
frequency lower than input frequency can be produced. 
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RF Cavity 



When an electron goes round the ring, it loses energy to 
radiation. When it passes through an rf cavity, it experiences 
an electric field in the forward direction. This gives it a push 
and replaces the energy that it has lost to radiation. 

Consider a metallic, rectangular box. Make a tiny hole on one 
side of the box, and stick a tiny wire through it. 
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RF Cavity 



Apply an oscillating voltage to this wire. Electric field from 
current would then induces charges on inner walls of box. 
There is a natural frequency of oscillation for field and charges. 
Field can resonate and grow stronger and stronger - if the 
voltage on that wire through the hole has the same frequency. 

This resonant frequency is obtained by solving the wave 
equation. For a rectangular box, the wave equation can solved 
using a product of sine or cosine functions. The result is a 
standing wave with frequency: 

c //ra7r\ 2 /n7T\ 2 . fl7T\ 2 

where l,m,n are integers. Choose one side, say 6, to be shorter 
than the other two. Choose m = l,n = 0,/ = 1. 
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RF Cavity 




This choice gives an oscillating electric field that is always 
parallel to the side b. It changes direction at the oscillation 
frequency, but must remain parallel to b. It is uniform along 
direction b. 
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RF Cavity 




Now make two small holes on the two opposite faces. Can pass 
an electron into one hole, and the field would accelerate the 
electron to the other hole! If AB is the desired path, just put 
the box in the path so that the path goes through the holes, 
and it will do the job of accelerating the particle. 

Note that a cavity that uses the wave this way is called standing 
wave cavity. There are other ways, e.g. travelling wave cavity. 
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Gaussian Units 



Gaussian units is based on cgs (centimetre-gram-second) units. 
It is commonly used in electromagnetism to make equations 
look simpler. But this also means that conversion to SI units is 
tricky because the equation must also be changed. 



c = 29,979,245,800 = 3 10 10 



Quantity 


Symbol 


31 unit 


Gaussian unit 


electric charge 


q 


1 C 


«(10" 1 c)Fr 


electric current 


i 


1 A 


<->(10- 1 c) Fr/s 


electric potential 


<p 


1 V 


«(10 8 c~ 1 )statV 


voltage 


V 


electric Held 


E 


1 V/m 


«(10 6 c" 1 )statV/cm 


magnetic induction 


B 


1 T 


«(10 4 )Gs 


magnetic field strength 


H 


1 A/m 


«H>(4TTlO- 3 }Oe 


magnetic dipole moment 


M 


1 Am 2 


<-» (10 3 ) erg/Gs 


magnetic flux 


*m 


1 Wb 


^(10 8 )Gscm 2 


resistance 


R 


1 0 


«-» (10 9 c~ 2 ) s/cm 


resistivity 


P 


1 Q m 


«(10 11 c~ 2 )s 


capacitance 


C 


1 F 


^(10~ 9 (^jcrn 


inductance 


L 


1 H 


<->(10 9 c"" 2 )s 2 /cm 



Name 


Gaussian units 


SI units 


Lorentz force 


F = q (e + iv x B^j 


F = q (E + v x B) 


Coulomb's law 


F = r 

r 2 


1 ?1<?2„ 

F — r 

4-Treo t 2 



Gauss's law 
(microscopic) 


V - E = Airp 


V • E = p/e 0 


Gauss's law for magnetism: 


V ■ B = 0 


V-B =0 


Maxwell-Faraday equation 
(Faraday's law of induction): 


VxE = "ft 


VxE -~at. 


Ampere-Maxwell equation 
(microscopic): 


„ „ 4tt IdE 
VxB=— J + - — 
c c ot 


„ „ 1 r3E 
VxB=^ 0 J + c2 ^ 



http: //en.wikipedia. org/wiki/Gaussian_units 



Google "Gaussian units" for details on conversions, 
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Synchrotron Radiation 



[In this topic, Gaussian units are used.] 

When an electron passes through a magnet, it bends and 
radiates energy. This is called synchrotron radiation. The rate 
of radiation loss for an electron when it revolves around a circle 
is 

P = -\7 4 a 2 [Gaussian units] 

(Gaussian units) where q is electron charge, a is acceleration 
and 7 is Lorentz factor. 

This is Lienard's formula. To derive this and understand the 
basic physics, first derive the Larmor radiation formula. 
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Synchrotron Radiation 




An electron is accelerated for a short time St to a velocity Sv. 
This "jerk" creates a kink in the electric field. This kink travels 
outward at the speed of light, and forms the radiation that is 
generated. 
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Synchrotron Radiation 




(Av)t 



To estimate this transverse electric field component, consider 
the kink at a time t. The kink is just a slight distortion of the 
electric field. The magnitude of the electric field along the kink 
remain approximately the same as the Coulomb field of the 
electron at rest: q/r 2 [Gaussian units]. The radial component 
of the field is also approximately the same, i.e. 

E r = \ [Gaussian units] 
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Synchrotron Radiation 



With the help of the triangle FGH enclosing the kink, and using 
similar triangles, we can see that 

E ± (Av)ts\nO 

— ^ = [Gaussian units] 

E r cAt 

Note that the distance AB is approximated as (/\v)t because 
the distance travelled during the initial At has been neglected. 
Substituting t = r/c, the transverse component can now be 
obtained: 

qas\nO . 
Ej_ = ~ — [Gaussian units] 

rc z 

where the acceleration a = Av/At. Using the Poynting vector 

S = — E x B [Gaussian units] 

and taking E = B (cgs units for plane wave, true if the particle 
is oscillating sinusoidally), so that 

S = —E 2 [Gaussian units] 
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Synchrotron Radiation 



the energy flux in direction 0 can be derived: 

1 q 2 a 2 s\n 2 0 
S = ~ o - [Gaussian units] 

Note that this result predicts that radiation is maximum in a 
direction perpendicular to the acceleration of the electron. The 
total power can now be obtained by integrating over all 
direction. This gives: 

o 2 2 

P = =— [Gaussian units] 

This is Larmor's equation. 



To derive the final expression, we show that when an electron 
goes round a circle at relativistic speed, the acceleration in 
Larmor's equation is increased by a factor of 7 2 . 

Suppose that in the lab frame, the electron is is moving in the 
x direction at a particular instance in time, and accelerating in 
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Synchrotron Radiation 



the y direction. In the rest frame, the electron velocity would 
be zero, but it would still be accelerating in the y direction. 
This is where Larmor's equation is valid. 

Then we do a Lorentz transformation to the lab frame. The 
simplest way to understand this is from the definition of 
acceleration: a = d 2 x/dt 2 . Because of time dilation, the dt 2 in 
the denominator gives the factor of 7 2 . 

Since the acceleration, a, gains a factor of the 7 2 when we 
transform to the lab frame, we can insert these into Larmor's 
equation and get 

P = -\7 4 a 2 [Gaussian units] 

This is Lienard's formula for the power from synchrotron 
radiation. 



46 



Vector Potential 



[We switch back to SI units here.] 



The static electric potential <j> is defined by 

E = -V0 

It is simplifies calculations and is interpreted as work done on a 
charge. The magnetic vector potential A is defined by 

B = V x A 

It also simplifies calculations but is not normally given a 
physical interpretation. It is often remembered as "the thing 
whose curl is magnetic field". It was historically interpreted as 
electromagnetic field momentum per unit charge. Lets see how. 



First, note that when the above expression is substituted into 
Faraday's law 

dB 

V x E = 

dt 



47 



Vector Potential 



we get 



^ ^ 9VxA ^ dA 

V x E = = V x 

dt \ dt 



This means that in time varying field, E has the additional term 
—dA/dt. So static electric potential must be generalised to 

dt 

To show that A is field momentum per unit charge, consider a 
long solenoid with a charge q outside. 
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Vector Potential 



/ 



\ 



/((((((»(((((((( 



qSE 




m 



SB 



If the magnetic field B changes over a short time, an electric 
field is generated around the solenoid according to Faraday's 
law: 

E.dl = -— [ B.dS 

C dt JS 

The vector potential definition can be written using Stoke's law 
as 

r A.dl = I B.dS 

'c Js 
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Vector Potential 



Comparing the two equations for a short time dt, 

Edt = -dA 

-d(qA) = qEdt = Fdt = d(rav) 

So an increase in particle momentum means a decrease in qA. 
(Just like an increase in particle energy means a decrease in 
potential energy <?</>.) 

So qA can be interpreted as field momentum. 

We know that total kinetic energy and potential energy is 
conserved. In the same way, we expect mv + gA to be 
conserved. So we can think of this as a generalised momentum 

p = mv + qA 
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Hamiltonian 



Hamiltonian method is useful for calculations on particle 
motion in accelerators. The Hamiltonian is sum of kinetic and 
potential energy 

H = T + V = -mv 2 + # = — + qcj) 

2 2m 

T must be expressed in terms of momentum p^. Since the 
system included both field and particle, pi should be the 
generalised momentum on the last slide 

p = rav + qA 

This suggests that 

(P-<?A) 2 

H = h q$ 

2m 

To check, must derive equation of motion from Hamilton's 
equations: 

dqi dH dpi dH 

dt dpi dt dqi 

51 



Hamiltonian 



In this case, pi is a component of p and xi is a component of x. 
Lets start with x components. First equation gives 

dx _ dH _ p x — qA x 

dt dp x rn 

mx = p — qA 

which is correct. Second equation gives: 

dp x _ d_ / (p - gA) 2 

dt dx \ 2m 

(p - qA) dA 

= —q . 

m dx 

Expanding 




d(mx + qA x ) (mx) dA dcj) 
— —q q — 

dt m dx dx 

dmv x dA x f.dA x ,dA v .dA z \ d<j> 

= —q q x h y \- z — q — 

dt dt \ dx dx dx J dx 
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Hamiltonian 



Using the full differential for A x 

dmv x ( . dA x . d A x . d A x ; d A x 

— —Q x 1- y 1- z 1 

dt \ dx dy dz dt 

(.dA x . .8A y . dA z \ 80 

— q x 1- y + z — q — 

V dx dx dx J dx 

= _ q ( dAx + <¥k\ + q ( dAz + <^A - qE - q 9Ax 
\ dy dx J V dx dz J 1 dt 

Using definitions B = V x A and E = -V<j> - dA/dt, 

dmv x 

dt — qv y B z - qv z B y + qE x 

dmv ^ ^ 
= gv X B + qE 

dt 

We have recovered the Lorentz force equation. This shows 
that the Hamiltonian is correct. 
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In Conclusion 

We have 

• reviewed the maths needed, 

• derived Maxwell's equations, 

• derived the wave equation, 

• derived waveguide modes, 

• derived synchrotron radiation, 

• reviewed vector potential, and 

• derived the Hamiltonian. 
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